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Abstract 

This paper describes advances in microwave frequency standards 

using laser-cooled atoms at BNM-SYRTE. First, recent improvements 

of the 133 Cs and 87 Rb atomic fountains are described. Thanks to 

the routine use of a cryogenic sapphire oscillator as an ultra-stable 

local frequency reference, a fountain frequency instability of 1.6 x 
10 -i4 r -i/2 where 

r is the measurement time in seconds is measured. 
The second advance is a powerful method to control the frequency 
shift due to cold collisions. These two advances lead to a frequency 
stability of 2 x 10~ 16 at 50 000 s for the first time for primary stan- 
dards. In addition, these clocks realize the SI second with an accuracy 
of 7 x 10~ 16 , one order of magnitude below that of uncooled devices. 
In a second part, we describe tests of possible variations of funda- 
mental constants using 87 Rb and 133 Cs fountains. Finally we give an 
update on the cold atom space clock PHARAO developed in collab- 
oration with CNES. This clock is one of the main instruments of the 
ACES/ESA mission which is scheduled to fly on board the Interna- 
tional Space Station in 2008, enabling a new generation of relativity 
tests. 
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1 Introduction: Einstein's legacy in modern 
clocks 

Modern clocks using laser cooled atoms owe a great deal to the famous 1905 
"annus mirabilis" of Einstein. Indeed the three theoretical problems that 
Einstein beautifully solved in 1905 are key ingredients in current atomic 
clocks, hundred years after Einstein's work. 

1. First the quanta of light, photons, are routinely used to cool atoms to 
microkelvin temperatures and to confine them in electromagnetic traps. 
Atom manipulation is a direct application of energy and momentum 
exchanges between light and matter. At one microkelvin, cesium atoms 
which form the basis for the current definition of the SI unit of time, the 
second, move at an average speed of 7mm.s _1 , enabling extremely long 
observation times and thus precision measurements. On Earth, atomic 
fountains enable unperturbed ballistic flight with duration approaching 
one second. Furthermore, every experiment in atomic physics routinely 
uses the Einstein's photoelectric effect in photodiodes to detect, and 
control light beams. Furthermore, the concept of photon is intimately 
connected to the famous Planck relationship E = hu between energy, 
Planck's constant and frequency of electromagnetic radiation, which is 
of paramount importance in atomic clocks. 

2. Second Einstein's theory of brownien motion with the famous relation- 
ship kftT = D/ot between temperature, diffusion coefficient and friction 
coefficient not only proved the existence of atoms, but beautifully ap- 
plies to Doppler and sub-Doppler laser cooling mechanisms at work in 
every cold atom experiment £Q . In optical molasses atoms are viscously 
confined by the bath of photons, they experience a three-dimensional 
random walk in position and storage times in excess of 10 seconds have 
been observed for this brownien motion. 

3. Third Einstein's theory of special (and later general) relativity intro- 
duced a new approach relating space and time, and the fundamental 
concept of relativistic invariance and Lorentz transformation. Einstein 
predicted that time in a fast moving frame seems to slow down to some- 
one not moving with it, and distances appear shorter. This revolution- 
ary approach had major fundamental as well as practical consequences 
in the following century. Clocks in different reference frames tick at dif- 
ferent rates and the well-known GPS receivers which equip boats, cars, 
and planes use routinely Einstein's relativity to determine their posi- 
tion with 10 meter accuracy. Indeed, the atomic clocks onboard the 24 
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GPS satellites orbiting the Earth at an altitude of 20 000 kms must be 
corrected for relativistic effects (time dilation and gravitational shift) 
in order to be synchronized with ground clocks and to reach this posi- 
tioning accuracy. The correction is about 38 microseconds per day. If 
each satellite would not apply this compensation the positioning error 
would reach 11 kilometers per day ! 

These three papers have had revolutionary consequences in Science and 
society. 

Historically, clocks have played a major role in tests of predictions of 
relativity theories, from the Hafele-Keating clock transport in jet planes, 
the Pound- Rebka gravitational shift measurement, the Vessot-Levine GP- 
A Space hydrogen maser red-shift measurement, and radar ranging Shapiro 
delay experiment [2|. In addition the current definition of time in the SI unit 
system, the second, relies on Einstein Equivalence Principle (EEP). This 
principle is the foundation for all gravitational metric theories that describe 
gravity as a consequence of curved space-time. The Einstein Equivalence 
Principle states 

1. if an uncharged test body is placed at an initial event in spacetime and 
given an initial velocity there, then its subsequent trajectory will be 
independent of its internal structure and composition 

2. In any freely falling frame, the outcome of any local non gravitational 
test experiments is independent of the velocity of the frame 

3. the outcome of any local non gravitational test experiment is indepen- 
dent of where and when in the universe it is performed 

An immediate consequence of EEP is that the fundamental constants of 
physics such as the gravitational constant G, or the fine structure constant 
a = e 2 /Aneohc must be independent of time and space. 

In this article we first describe recent progress in the realization of the SI 
second using laser cooled cesium and rubidium clocks. In a second part we 
use these highly stable devices to perform new tests of Einstein Equivalence 
Principle, namely the constancy of fundamental constants. 

2 Atomic Fountains 

The ever increasing control of the motion of atomic samples is at the origin 
of recent progress in atomic frequency standards and precision measurements 
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[3]. Laser cooled and trapped atoms enable long observation times required 
for high precision measurements. Charged particles confined in Paul or Pen- 
ning traps offer extremely long storage enabling high precision mass mea- 
surements, fundamental tests, and the realization of ultra-stable microwave 
and optical clocks. The recent NPL frequency measurement of an optical 
transition in Sr + ion with an uncertainty of 3 10~ 15 [3] is only a factor three 
or four worse than the current accuracy of cesium fountains. Precision mea- 
surements with neutral atoms on the other hand are usually performed in 
an atomic fountain where laser cooled atoms ballistically propagate for dura- 
tions up to one second. In the last decade, atomic clocks and inertial sensors 
using matter wave interferometry in fountains have become two of the most 
important applications of cold atoms [3J E] . About two dozens of fountain 
devices are now used for a variety of applications. It has been shown recently 
that microwave and optical clocks as well as matter-wave inertial sensors be- 
long to the same general class of atom interferometers [Hj . As an example the 
current sensitivity in acceleration measurement with atom interferometers is 
on the order of 3 x 10~ 8 m.s -2 in one minute measurement duration and, in a 
decade, cesium fountain clocks have gained almost two orders of magnitude 
in accuracy. As we show in this paper the fractional inaccuracy of the BNM- 
SYRTE fountains at Paris Observatory do not exceed today 7 x 10~ 16 which 
corresponds to less than a single second error over 50 million years, allowing 
for the realization of SI unit of time, the second, at the same level. About 
half a dozen fountains throughout the world at metrology institutes including 
PTB, NIST, IEN, NPL, have now an accuracy near 1CT 15 , making fountains 
a major contributor to the accuracy of the TAI (Temps Atomique Interna- 
tional). In the future, many applications, such as positioning systems (GPS, 
GALILEO, GLONASS) as well as scientific applications will benefit from 
these developments. For instance, deep space satellites have travel durations 
of several years across the solar system. Precise monitoring of their position 
requires timescales with very low long term drift. Also, using advanced time 
and frequency transfer systems (operating at higher carrier frequency and 
chip rate and/or using two way transfer techniques) may lead to positioning 
accuracy at the millimeter level for averaging time of a few hundred seconds. 
This would impact many geodetic applications. 

In this paper we show that prospects for further improvements are im- 
portant. A frequency comparison between two fountains exhibits a stability 
of 2 10 -16 at 50 000 second averaging time, for the first time for atomic stan- 
dards. This frequency resolution sets the stage for clock accuracy at the 
10~ 16 level for cesium, almost one order of magnitude potential gain, and 
even better for rubidium with its far reduced collision shift IE] . We be- 
gin by recalling the basic operation of fountain atomic clocks and introduce 
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several new techniques which demonstrate frequency measurements with a 
frequency resolution at the 10 -16 level. The first technique makes use of 
an ultra-stable cryogenic oscillator to interrogate the clock transition in the 
fountain. Thanks to its extremely good short term frequency stability and 
low phase noise, the frequency stability of cesium and rubidium fountains is 
one order of magnitude below that of fountains using an ultra-stable quartz 
oscillator as interrogation oscillator. It currently reaches 1.6 x 10~ 14 r -1 / 2 
where r is the averaging time in seconds. The fundamental quantum noise 
of the clock is now reached with atomic samples of up to 10 7 atoms. The 
second advance deals with a new technique to measure and cancel with high 
precision the collisional shift in the clock. This shift is a major plague in 
cesium clocks and is much reduced (two orders of magnitude) in rubidium 
devices |7J |H]. The method uses interrupted adiabatic population transfer 
to prepare precise ratios of atomic densities. We show here that the cesium 
collisional shift can be measured and cancelled at the 10~ 3 level. By com- 
paring rubidium and cesium fountains over a duration of six years, a new 
upper limit for the drift of fundamental constants has been obtained. Fi- 
nally we present the development status of the PHARAO cold atom space 
clock which is under industrial realization. PHARAO will fly onboard the 
international Space station in the frame of the European ACES mission in 
2008-2009 and perform fundamental physics tests such as an improved mea- 
surement of Einstein's red-shift, search for drift of fundamental constants 
and special relativity tests. 



3 Recent advances in cesium and rubidium 
fountains 

In this section we briefly review recent advances on cesium and rubidium 
fountains performed in our laboratory, BNM-SYRTE where three laser cooled 
atomic fountains are in operation. The first one (FOl), in operation since 
1994 !5J, has been refurbished recently. The second one (FOM), a trans- 
portable fountain, is derived from the PHARAO space clock prototype p. 
This fountain was transported on two occasions to the Max Planck Insti- 
tute in Garching for direct frequency measurement of the hydrogen Is — > 2s 
transition |U|(See section on stability of fundamental constants). The third 
one (F02), a dual fountain operating with 133 Cs or 87 Rb, is described in [7]. 
Here we only briefly describe the present design and recent improvements of 
FOl and F02. A scheme of the fountain apparatus is shown in Fig[Q An 
optical bench provides through optical fibers all beams required for manipu- 
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lating and detecting the atoms. The fountains operate with lin _L lin optical 
molasses. Atoms are cooled by six laser beams supplied by preadjusted fiber 
couplers precisely fixed to the vacuum tank and aligned along the axes of 
a 3 dimensional coordinate system, where the (111) direction is vertical. In 
FOM, optical molasses is loaded from a 133 Cs vapor and 3 x 10 7 atoms are 
cooled in 400 ms. In FOl and F02, optical molasses are loaded from a 
laser slowed atomic beam which is created by diffusing 133 Cs or 87 Rb vapor 
through a bundle of capillary tubes. With this setup 3 x 10 8 133 Cs atoms can 
be loaded in 400 ms in FOl. In F02 an additional transverse cooling of the 
atomic beam increases the loading rate to 10 9 atoms in 100 ms for 133 Cs. 
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Figure 1: Schematic view of the dual Cs-Rb atomic fountain. 

The atoms are launched upwards at 4 m.s -1 by using moving optical mo- 
lasses and cooled to ~ 1 /iK in the moving frame by adiabatically decreasing 
the laser intensity and increasing the laser detuning. In normal operation 
atoms in the clock level \F — 3, rriF = 0) are selected by microwave and light 
pulses. 
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Figure 2: BNM-SYRTE fountain ensemble. 

About 50 cm above the capture zone, a cylindrical copper cavity (TE n 
mode) is used to probe the hyperfme transition in a Ramsey interrogation 
scheme. The cavities have a loaded quality factor of Qfoi = 10000 for FOl 
and Qfo2 — 6600 for F02. Both cavities can be fed with two coupling 
irises oppositely located on the cavity diameter. Symmetric or asymmetric 
feedings are used to evaluate and reduce the residual Doppler effect due to 
imperfections of the standing wave in the cavity and a tilt of the launch 
direction of the atoms. 

The microwaves feeding the cavities are synthesized from the signal of an 
ultra-stable cryogenic sapphire resonator oscillator (CSO) developed at the 
University of Western Australia As shown in FigEl the three fountains 
use the same CSO oscillator to synthesize the microwave signals probing 
the atomic transition. To reduce its drift, the CSO is weakly phase-locked 
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Figure 3: Experimental Ramsey fringes (transition probability as a function 
of the microwave detuning) measured with 133 Cs in F02 fountain. The insert 
shows the central fringe with a FWHM of ~ 1 Hz. Each point is a single 1.3 s 
measurement. At half maximum of the central fringe, the signal to noise ratio 
is 5 000, within 20% of the fundamental quantum noise with ~ 10 7 detected 
atoms. 

to a hydrogen maser. This maser contributes to the local timescale and 
TAI (Temps Atomique International) through various time and frequency 
transfer systems. With this setup, atomic fountains are used as primary 
frequency standard to calibrate TAI and can be compared to other remote 
clocks. Nowadays, atomic fountains are the dominant contributors to the 
accuracy of TAI. 

The 11.932 GHz output signal from the CSO is converted in order to 
synthesize 11.98 GHz and 100 MHz signals, both phase coherent with the 
H-maser. F02 uses the 11.98 GHz signal to generate 9.192 GHz by a home- 
build low noise synthesizer which achieves a frequency stability of 3 x 10~ 15 
at 1 s by operating only in the microwave domain. This scheme reduces 
at the minimum the phase noise and the spurious side-bands induced by the 
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down conversion process. A similar setup is used to synthesize the 6.834 GHz 
required for the F02 fountain operation with 87 Rb. The 150 m distance be- 
tween FOl, FOM and the CSO prevents the direct use of the 11.98 GHz sig- 
nal. Instead, a 100 MHz signal is synthesized from the CSO and distributed 
to FOl and FOM via a high stability RF cable. Finally, a 100 MHz to 
9.192 GHz home-made synthesizer generates the interrogation signal. These 
additional steps degrade the phase noise of the interrogation signal in FOl 
and FOM with a frequency stability currently limited to ~ 2 x 10~ 14 at 1 s. 

Frequency stability 



tt 1 1 — i i i i i 1 1 1 1 — i i i i 1 1 1 1 1 — i — i i i 1 1 1 1 1 — i — i i i 1 1 1 1 r 




time (s) 

Figure 4: Fractional frequency instability of F02 against CSO for high den- 
sity (HD, red squares) and low density (LD, green circles) configurations. It 
demonstrates a stability of 1.6 x 10~ 14 r~ 1//2 for a 133 Cs primary standard. 
Also shown is the fractional frequency instability for the differential mea- 
surement between HD and LD (blue triangles). This curve demonstrates an 
excellent rejection of the CSO fluctuations in the differential measurement, 
allowing for a fractional frequency resolution of 2.5 x 10~ 16 at 20 000 s. In 
this measurement, the collisional shift at LD is the frequency difference be- 
tween HD and LD (~5 x 10~ 14 ). It is obtained with a resolution close to 2 
parts in 10~ 16 and it is stable at the 0.5% level over 20 000 s. 

Atoms selected in \F = 3, mp = 0) cross the microwave cavity on the way 
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up and on the way down, completing the two Ramsey interactions. After the 
Ramsey interrogation, the populations N e and N g of both clock levels |e) 
and \g) are measured by fluorescence. The number of detected atoms is typ- 
ically 0.5% of the initially captured atoms. The signal p = N e /(N e + N g ) is 
equal to the atomic transition probability and is insensitive to atom number 
fluctuations. A typical Ramsey resonance is presented in FigOJl From the 
transition probability, measured on both sides of the central Ramsey fringe, 
we compute an error signal to lock the microwave interrogation frequency to 
the atomic transition using a digital servo loop. At the quantum limit one 
expects S/N = l/os p = 2yN for iV detected atoms, where a$ p is the shot 
to shot standard deviation of the fluctuations of the transition probability. 
The frequency corrections are applied to a computer controlled high resolu- 
tion DDS synthesizer in the microwave generator. These corrections are used 
for the accuracy and frequency stability evaluations of each fountain. The 
fractional frequency instability of the F02 fountain operating with ~ 10 7 de- 
tected atoms and measured against the cryogenic oscillator is plotted FigEJ 
At the quantum limit one expects a frequency instability, characterized by 
the fractional Allan standard deviation, given by: Oy(r) = (l/TcQ at )^T c /NT, 
where Q at ~ 10 10 is the atomic quality factor, r and T c are respectively the 
averaging time and the cycle duration. Above the servo-loop time constant 
(~ 3 s) and below 100 s, the fractional instabilities of FOl and F02 are 
(7 y (r) = 2.9 x lO-^r" 1 / 2 and 1.6 x lO" 14 ^ 1 / 2 respectively, within ~ 20% 
of the standard quantum limit. For longer averaging time the frequency 
instability is dominated by the frequency fluctuations of the CSO and the 
H-maser. This is the first demonstration of routinely operated primary fre- 
quency standards with frequency instabilities in the low 10 _14 r _1//2 region. 
We will show below that this excellent short term stability enables an evalua- 
tion of systematic frequency shifts and frequency comparisons between clocks 
at the 10~ 16 level in a few days. 



Accuracy 

All known systematic frequency shifts are evaluated in our fountains. The 
accuracy budget for each shift is given in table I for 133 Cs. The overall 
uncertainty, the quadratic sum of all uncertainties is 7.5 x 10~ 16 for FOl, 
6.5 x 10~ 16 for F02 and 8 x 10~ 16 for FOM. In the following, we only discuss 
some of the most bothersome effects and the recent improvements in their 
evaluation. A more complete discussion of systematic effects can be found in 

nn. 
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Table 1: Systematic fractional frequency shifts for FOl and F02. 





FOf (xfO 16 ) 


F02 (xlO 16 ) 


FOM (xlO 16 ) 


Quadratic Zeeman effect 


1199.7 ±4.5 


1927.3 ±0.3 


351.9 ±2.4 


Blackbody radiation 


-162.8 ±2.5 


-168.2 ±2.5 


-191.0 ±2.5 


Collisions and cavity pulling (HD) 


-197.9 ±2.4 


-357.5 ±2.0 


-34.0 ±5.8 


Spectral purity & leakage 


0.0 ±3.3 


0.0 ±4.3 


0.0 ±2.4 


First order Doppler effect 


< 3 


< 3 


< 2 


Ramsey & Rabi pulling 


< 1 


< 1 


< 1 


Microwave recoil 


< 1.4 


< 1.4 


< 1.4 


Second order Doppler effect 


< 0.08 


< 0.08 


< 0.08 


Background collisions 


< 1 


< 1 


< 1 


Total uncertainty 


±7.5 


±6.5 


±7.7 



Cold collisions and cavity pulling 

The cold collision frequency shift is known to be particularly large for 133 Cs 
O HB] . For instance, when F02 is operated at its best frequency stability 
the shift is ~ 10~ 13 . The linear extrapolation of this effect to zero density is 
known to be delicate. As pointed out in the first paper observing the cold 
collision shift in 133 Cs fountains ^2], selecting atoms in the clock levels using 
microwaves may lead to distortions of the position or velocity distribution. 
Methods to cope with these effects have been proposed [SJ, yet the linear 
extrapolations have proved to be valid only at the 5% to 10% level. 

To evaluate the collision shift at the 10~ 3 level (a requirement for a fre- 
quency stability and accuracy at 10~ 16 ), we recently developed a new method 
based on interrupted adiabatic passage to select atoms in the \F = 3, m-p = 0) 
state [H]. Atomic samples are prepared by transferring atoms from \F = 
4,m F = 0) to \F = 3,m F = 0) with an efficiency precisely equal to 100% 
(high density, HD) or 50% (low density, LD). With this method, the atom 
number is changed without affecting either the velocity or the position dis- 
tributions. Therefore, the density ratio LD/HD is equal to the atom number 
ratio and is 1/2 at the 10 -3 level. Since the collisional shift is proportional to 
the atomic density, it can be extrapolated to zero density with this accuracy. 
In addition, with this method, the cavity frequency pulling |Hl EH] is also 
accounted for. 
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Number of cycles 



Figure 5: Fractional instability of the ratio of the detected atom number in 
\F = 4,rriF = 0) between low density and high density configurations as a 
function of the number of fountain cycles. The measured ratio is 0.5005(2). 
Each cycle lasts ~ 1.3 s. The stability (solid line) decreases as the square 
root of the number of cycles. 

The collisional shift is measured in real-time with the following differential 
method. The clock is operated alternately in the HD configuration for 60 s 
and in the LD configuration for the next 60 s. This timing choice minimizes 
the noise due to frequency instabilities of the CSO oscillator. As seen in 
FigEJ at 120 s the stability of F02 against CSO is near its minimum. Also, 
over this 120 s period, density fluctuations do not exceed ~ 1%. One the 
other hand, due to slow changes in the clock environment, we observe that 
the density may fluctuate up to 10 — 20% over one or several days. Our 
differential method efficiently cancels these slow daily density variations. 

In [Hj, our calculations predicted that the interrupted adiabatic passage 
method does provide a LD/HD ratio precisely equal to 1/2 to better than 
10 -3 . Initially, we were experimentally able to realize this ratio at the 1% 
level. Improvements in the accuracy of the microwave frequency synthesis 
for the adiabatic passage enable us to now reach a precision of 2 x 10~ 3 for 
this ratio. 

During routine operation of the fountains, the number of detected atoms 
in each hyperfine state is recorded for both LD and HD configurations. As 
seen in Fig|3 the Allan standard deviation of the measured LD/HD atom 
number ratio decreases as the square root of the number of fountain cycles 
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(or time), down to a few parts in 10 4 for one day of averaging. Despite the 
10 — 20% slow drift in atom number over days, this ratio remains remarkably 
constant. The LD/HD atom number ratio in \F — 4, mp = 0) is found equal 
to 1/2 to better than 10~ 3 . This method relies on fluorescence measurements 
made in the detection zones for each fountain cycle. Various measurements 
have been performed to establish that the measurement of this ratio is not 
biased by more than 10~ 3 due optical thickness effects in the detection. On 
the other hand, the LD/HD atom number ratio in \F = 3, mp = 0) is found 
to slightly differ from 1/2 by 0.3% typically. This deviation originates from 
atoms in the \F = 3,mF ^ 0) states populated by imperfections in the state 
preparation. This deviation must be taken into account in the evaluation 
of the collisional shift. In [HQ, we showed that the frequency shift of the 
clock transition due to \F = 3, m-p ^ 0) atoms is at most 1/3 of that of 
collisions between \F = 3,mp = 0) and \F = 4, mp = 0) clock states. Their 
contribution to the collisional frequency shift is thus at the 0.1% level. In 
summary, when this adiabatic passage method is used, we take a 2 x 10 -3 
relative uncertainty for the determination of the high density cold collision 
shift. 

Effect of microwave spectral purity and leakage 

Spectral impurities of the interrogation signal and microwave leakage may 
cause shifts of the clock frequency. In order to evaluate these effects, we 
make use of their dependance with the microwave power. We alternate the 
Ramsey interrogation between a configuration of tt/2 and 37r/2 pulses, i.e. 
a variation of a factor of 9 in microwave power. Within the resolution of 
the measurement of 3.3 x 10 -16 , no frequency shift is observed. In this 
measurement, four data sets are recorded, LD and HD at tt/2 and LD and 
HD at 37r/2. In this way, the collisional shift (which may also change with 
the microwave power) is evaluated and cancelled for both tt/2 and 37r/2 
configurations by the differential method described above, allowing for the 
extraction of a possible influence of microwave spectral purity and leakage 
alone. 

Residual first order Doppler effect 

A frequency shift due to the first order Doppler effect can occur if the mi- 
crowave field inside the interrogation cavity exhibits a phase gradient and 
the atoms pass the cavity with a slight inclination from the cavity axis. We 
determine the frequency shift due to the linear component of the phase gra- 
dient in a differential measurement by coupling the microwave interrogation 
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signal "from the left", "from the right" or symmetrically into the cavity, 
providing 3 data sets. The observed shift between the "left" and symmetric 
configuration is (—25.3 ±1.1) x 1CT 16 while the shift between the "right" 
and symmetric configuration is (+24.0 ± 1.2) x 10 -16 . The magnitude of 
this residual first order Doppler effect is consistent with a simple estimate 
of the residual traveling wave component in the cavity J7j together with a 
misalignment between the local gravity and the launch direction < 1 mrad. 
The mean of these two measurements is (—0.7 ± 0.8) x 10~ 16 and consistent 
with zero, indicating that the traveling wave component is well cancelled 
in the symmetric coupling configuration. Using the atoms as a probe, we 
can indeed ensure that the cavity is fed symmetrically to better than 1% 
in amplitude and 60 mrad in phase, which cancels the effect of linear phase 
gradient to ~ 1%, better than the above measurement resolution. As a conse- 
quence, only the quadratic phase dependence of the microwave field remains 
as a possible source of residual Doppler shift. A worst case estimate based 
on P?] gives an upper bound for the fractional frequency shift of 3 parts in 
10 16 , which we conservatively take as the overall uncertainty associated with 
residual first order Doppler effect. 

Other contributions to the accuracy budget are listed in Table The 
total accuracy currently reaches 7.5 parts in 10 16 for FOl, 6.5 parts in 10 16 
for F02, and 8.0 parts in 10 16 for FOM. This represents a one order of 
magnitude improvement over uncooled cesium devices. In the future, we 
anticipate that the extensive use of the methods described above will enable 
us to bring the accuracy of 133 Cs fountains below 2 parts in 10 16 and the 
accuracy of 87 Rb to a even lower value thanks to its 100-fold lower collision 
shift 0|H!. 

Frequency comparisons between two 133 Cs fountains at 2 10 16 

The routine operation of two atomic fountains near the quantum noise limit 
using the CSO as an interrogation oscillator enables frequency comparisons 
in the low 10~ 16 range, for the first time for primary frequency standards. 
Figure |H1 presents the frequency stability between FOl, F02 and CSO. Each 
fountain is operated in differential mode in order to permanently evaluate 
and cancel the collision shift. Appropriate post-processing of the data thus 
enables us to construct for each fountain, a clock which is free of the cold 
collision shift and whose stability is shown in Fig. |H1 against the CSO oscilla- 
tor. Fig. El also shows that the combined stability between these two clocks 
reaches 2.2 x 10~ 16 at 50 000 s, a previously unattained long term stability. 
From this data, we infer that at least one of the two fountains has a stability 
below {2.2/ y/2) x 10~ 16 = 1.6 x 10" 16 at the same averaging time. The mean 
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Figure 6: Fractional frequency instability (Allan deviation) between FOl 
and F02 fountains (blue triangles). After 50 000 s of averaging, the stability 
between the two fountains is 2.2 parts in 10 16 . Also plotted is the fractional 
frequency instability of FOl (red circles) and F02 (black squares) against 
the CSO locked to the hydrogen maser. 

fractional frequency difference between the two fountains is 4 x 10~ 16 , fully 
compatible with the accuracy of each of the two clocks as stated in Table H 
This very good stability sets a new challenge for time and frequency trans- 
fer systems between remote clocks. As an example, long distance frequency 
comparisons between PTB and NIST fountains were performed at the level 
of only 6 x 10 -16 after 2 weeks of averaging with GPS jTH]- Similarly, com- 
parisons between BNM-SYRTE and PTB recently achieved 2 x 10~ 15 for one 
day of integration with TWSTFT P]. 

4 Einstein Equivalence Principle and Stabil- 
ity of fundamental constants 

Highly accurate atomic clocks offer the possibility to perform laboratory 
tests of possible variations of fundamental constants. Such tests interestingly 
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complement experimental tests of the Local Lorentz Invariance and of the 
Universality of free-fall to experimentally establish the validity of Einstein's 
Equivalence Principle (EEP). They also complement tests of the variability 
of fundamental constants on different timescales, geological timescale [201 I2J| 
and cosmological timescale [221 EHj Nearly all unification theories (in par- 
ticular string theories) violate EEP at some level |2"H |23 |2H] which strongly 
motivates experimental searches for such violations. 

Tests described here are based on highly accurate comparisons of atomic 
energies. In principle, it is possible to express any atomic energy as a function 
of the elementary particle properties and the coupling constants of funda- 
mental interactions using Quantum Electro-Dynamics (QED) and Quantum 
Chromo-Dynamics (QCD). As a consequence, it is possible to deduce a con- 
straint to the variation of fundamental constants from a measurement of the 
stability of the ratio between various atomic frequencies. 

Different types of atomic transitions are linked to different fundamental 
constants. The hyperfine frequency in a given electronic state of alkali-like 
atoms (involved for instance in 133 Cs, 87 Rb [23, 199 Hg+ |2H1 [21, 171 Yb+ 
microwave clocks) can be approximated by: 

l^cxi^/)^] « 2 F«(«), (1) 

where the superscript (i) indicates that the quantity depends on each partic- 
ular atom. Roo is the Rydberg constant, c the speed of light, grW the nuclear 
g-factor, m e /m p the electron to proton mass ratio and a the fine structure 
constant. In this equation, the dimension is given by -RqoC, the atomic unit 
of frequency. A^ s is a numerical factor which depends on each particular 
atom. F^ s (a) is a relativistic correction factor to the motion of the valence 
electron in the vicinity of the nucleus. This factor strongly depends on the 
atomic number Z and has a major contribution for heavy nuclei. Similarly, 
the frequency of an electronic transition (involved in H j^HI, 40 Ca [ST], 199 Hg + 
|32*| . 171 Yb + [SSI Ell optical clocks) can be approximated by 

i/t ~ R^c x x F»>). (2) 

Again, the dimension is given by RooC. A^ ec is a numerical factor. F^ c (a) 
is a function of a which accounts for relativistic effects, spin-orbit couplings 
and many-body effects 1 . 

1 It should be noted that in general the energy of an electronic transition has in fact a 
contribution from the hyperfine interaction. However, this contribution is a small fraction 
of the total transition energy and thus carries no significant sensitivity to a variation of 
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According to 



the sensitivity to g-factors g^> and to the proton 



mass m p can be related to a sensitivity to fundamental parameters, namely 
the mass scale of QCD Aqcd and the quark masses m q = {m u + m<z)/2 and 
m s . Therefore, any measurement of the ratio between atomic frequencies can 
be interpreted as testing the stability of 4 dimensionless fundamental con- 
stants: a, m 9 /A QCD , m s /A QCD and m e /A QCD . The sensitivity to m s /A QCD 
is relatively weak compared to the 3 other constants. The sensitivity coeffi- 
cients have now been calculated for a large number of atomic species used in 
atomic clocks ESI EZ1 EHl ISHl E01 EH EEl E2J- Reliable knowledge of these 
sensitivity coefficient at the 1% to 10% level is required to deduce limits to 
a possible variation of each of these fundamental parameters by combining 
the results of several complementary clock comparisons. 
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Figure 7: Measured 87 Rb frequencies referenced to the 133 Cs foun- 
tains over 72 months. The 1999 measurement value (i/Rb(1999) = 
6 834 682 610.904 333 Hz) is conventionally used as reference. A weighted lin- 
ear fit to the data (solid line) gives f t In = (-0.5 ± 5.3) x lO^yr" 1 . 
MJD represent Modified Julian Dates. 

Figd summarizes the comparison between 87 Rb and 133 Cs hyperfine fre- 

fundamental constants. The same applies to higher order terms in the expression of the 
hyperfine energy JQ). A precision of 1 to 10 % on the sensitivity is sufficient to interpret 
current experiments. 
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quencies that have been performed using the above described fountain en- 
semble over a duration of 6 years. Each point on the graph summarizes the 
result of one to two months of measurements which include each time an 
evaluation of all known systematic effects fZT \ I43 [ 14*4*] . A weighted linear fit 
to the data in Fig0 determines how these measurements constrain a possible 
time variation of ^Rb/t'cs- We find: 

* in = (-0.5 ± 5.3) x 1(T 16 yr" 1 (3) 

which represents a 100-fold improvement over the Hg + -H hyperfine energy 
comparison [28 . This results implies the following constraint: 

" i_ ( 9Cs n.49 \ _ /n c _u c o\ ^ in-16,^-1 



In £h± a o.4» = (o.5 ± 5.3) x 10" lb yr" 1 . (4) 
dt \g^h J 

Using the calculated link between g-factors and m q , m s and Aqcd (ref. 
jHSl EE]), we find the following constraint to the variation of fundamental 
constants: 

j t In (a ' 49 K/A Q cD] am K/Aqcd] - 027 ) = (0.5 ± 5.3) x lO" 16 yr" 1 . (5) 

As pointed out in [2SJ E3 EE] , the hypothetical unification of all interactions 
implies that a variation of the fine-structure constant a should be accompa- 
nied by a variation of the strong interaction constant and of elementary par- 
ticle masses. Within this framework, current estimates gives |23] EE1 EH HE] : 

S(m/A Q cp) 5a 

- — — ~ 35 x — . (6) 

(m/ Aqcd) a 

Within this theoretical framework, the present comparison between Rb and 
Cs fountains (Eq. E]) constrains a time variation of a at the level of 7 x 
10 -1 ' yr -1 . In the future, improvement of 87 Rb and 133 Cs fountain to ac- 
curacies of few parts in 10 16 and repeated comparisons over several years 
between these devices will improve the above result by at least one order of 
magnitude. 

The transportable fountain FOM has similarly been used as a primary 
standard in the measurement of the frequency of the hydrogen 1S-2S 
transition performed at Max-Planck Institut in Garching (Germany) [3lH l47j. 
Two measurements performed over a 4 year period constrain fractional vari- 
ations of t'cs/^H at the level of (3.2 ± 6.3) x 10~ 15 .yr -1 . This constrains 
fractional variations of gcs(^e/ 'm p )a 2 ' 83 at the same level [2H1ISZ]- Combin- 
ing these results with other recent comparisons ( 199 Hg + optical clock versus 
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133 Cs fountain OHl, 171 Yb+ optical clock versus 133 Cs fountain 
it is possible to independently set limits on variations of a, gRb/gcs and 
gcs{m e /m p ). These measurements test the stability of the electroweak in- 
teraction (a) and of the strong interaction (gnb/gcs, PcsO^eA^p)) separately 
(13 EH] and independently of any cosmological model. 

5 The PHARAO cold atom space clock and 
ACES 

PHARAO, as "Projet d'Horloge Atomique par Refroidissement d'Atomes en 
Orbite", has started in 1993 with the objective of performing fundamental 
metrology with a space cold atom clock jHOj- The combination of laser cooling 
techniques [51 J and microgravity environment indeed allows the development 
of space clocks with unprecedented performances. 

To demonstrate the feasibility of a compact cold atom clock operating 
in microgravity, BNM-SYRTE and LKB with the support of CNES (the 
French space agency) have undertaken the construction of a clock prototype 
in 1994. The prototype was successfully tested in 1997 in jet plane parabolic 
flights [H]. The same year, ESA, the European Space Agency, selected the 
ACES proposal (Atomic Clock Ensemble in Space) [52]. ACES will perform 
fundamental physics tests by using the PHARAO cold atom clock, a H-maser 
(developed by the Neuchatel Observatory) and a time and frequency transfer 
system MWL on a platform developed by ESA. This ensemble will fly on 
board the international space station in 2008-2009. The station is orbiting 
at a mean elevation of 400 km with a 90 mn period and an inclination angle 
of 51.6°. The planned mission duration is 18 months. During the first 6 
months, the performances of the PHARAO cold atom clock in space will be 
established. Thanks to the microgravity environment the linewidth of the 
atomic resonance will be varied by two orders of magnitude (from 11 Hz 
to 110 mHz). The target performance is 7 x 10 _14 r -1 / 2 for the frequency 
stability and 10 -16 for the frequency accuracy. In the second part of the 
mission, the onboard clocks will be compared to a number of ground based 
clocks operating both in the microwave and the optical domain. 

In 2001, PHARAO entered into industrial development under the man- 
agement of CNES with the construction of two clock models, an engineering 
model for test and validation, and a flight model. 
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Figure 8: The PHARAO sub-systems and interfaces. 
The PHARAO instrument 

The clock is composed of four main sub-systems as shown in figure |H1 Each 
sub-system has been subcontracted to different manufacturers and they will 
be assembled at CNES Toulouse to validate the clock operation. 

The laser source provides all the laser tools for cooling, launching and 
detection of the atoms. Two extended cavity diode lasers [SB] are used as 
master lasers. One of them injection-locks two slave diode lasers to provide 
high laser power for capturing ~ 10 8 in optical molasses. The second laser is 
used as a repumping laser. The two master laser frequencies are stabilized by 
servo-loops using absorption signals through cesium cells. The other laser fre- 
quencies are synthesized by using 6 acousto-optic modulators (AOM). These 
AOM also control the laser beam amplitudes. The output laser beams are 
connected to the cesium tube through polarization maintaining optical fibers. 
Figure El shows the PHARAO optical bench during the assembly. The total 
mass is 20 kg, the volume is 26 liters and the power consumption is 40 W. 

The cesium tube provides the atomic source, the controlled environment 
for the atomic manipulation, the interrogation and detection process (figure 
EJ). Its design is similar to atomic fountains except for the interrogation zone 
where a two zone Ramsey cavity is used. The Ramsey cavity (figure ITT]) has 
been specially developed for this application and forms a ring resonator. One 
coupling system feeds two symmetrical lateral waveguides which meet at the 
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Figure 9: The PHARAO optical bench during assembly. The bench surface 
is 55 x 33 cm 2 . The laser source includes 8 frequency stabilized diode lasers. 
The laser light for atom manipulation is coupled to the vacuum tube through 
optical fibers, (courtesy of EADS SODERN). 

two interaction zones. The advantage of this configuration is to provide very 
weak phase disturbances of the internal microwave field while enabling large 
holes (8x9 mm) for the atom path. The flight model of the microwave cavity 
is currently mounted (September 2004) inside the atomic fountain FOl to 
measure the end to end cavity phase shift before integration in the flight 
model. These measurements and numerical simulations, should enable us to 
determine the cavity phase shift effect with an accuracy of a few parts in 
10 17 . 

The cesium tube is designed for a vacuum of 10~ 8 Pa in order to minimize 
the cold atom losses with the background gas collisions. Three layers of 
magnetic shields and a servo system maintain the magnetic field instability in 
the interaction zone below 20 pT. Similarly, the interaction zone temperature 
is regulated to better than 0.2°C. 

The microwave chain synthesizes the two radiofrequency signals for the 
state selection cavity and the interrogation cavity. A 100 MHz VCXO (Volt- 
age Control Oscillator) is phase-locked to an Ultra Stable Oscillator (USO) 
for the short term stability and to the Space Hydrogen Maser (SHM) for the 
medium term. Three USOs have been space qualified for our application. 
We have compared these quartz oscillators to the BNM-SYRTE CSO. Their 
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Figure 10: Cross-section of the cesium tube. The mass is 45 kg and the 
volume 70 liters. 

frequency stability is on the order of 7 x 10~ 14 from 1 to 10 s integration 
time. The engineering model of the chain has been fully tested and the re- 
sults are in agreement with the performance objectives of the space clock. 
A further performance verification is currently being made by using the mi- 
crowave source in the F02 fountain. All PHARAO sub-systems are driven 
and controlled by a computer (UGB, On Board Data Processing Unit). The 
UGB also manages the data flux between the clock and the ACES payload. 
When assembled, the clock fills a volume of about 200 1 for a weight of 91 kg 
and an electric consumption of 114 W. 

The final assembly of the engineering model of the PHARAO clock will 
start at the end of 2005 at CNES-Toulouse. After the clock functional and 
performance test are made, the flight model will be assembled and finally 
tested. For both models, we expect to reach 10~ 15 frequency accuracy in the 
Earth gravity environment and 10~ 16 in microgravity environment. 

Scientific objectives of the ACES mission 

The objectives of PHARAO/ACES are (i) to explore and demonstrate the 
high performances of the cold atom space clock (ii) to achieve time and 
frequency transfer with stability better than 10~ 16 and (iii) to perform fun- 
damental physics tests. A detailed account can be found in [El] . 

The combination of PHARAO with SHM will define an on board fre- 
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Figure 11: The Ramsey interrogation cavity. The cavity in pure copper is 
screwed on a rigid structure to avoid deformation during launch. The length 
of the cavity is 280 mm. The atoms enter the cavity through the cut-off 
waveguide with a rectangular shape (on the left). Also visible in the center 
is the the microwave coupling system. (Courtesy of EADS SODERN and 
TAS). 

quency reference having a long term stability and accuracy provided by 
PHARAO and a short term stability determined by SHM. The resulting 
fluctuations of ACES frequency reference are expected to be about 10 ps 
per day. The orbit of ISS will allow ground users to compare and syn- 
chronize their own clock to ACES clocks, leading to a worldwide access to 
the ultra stable frequency reference of ACES. The results of these compar- 
isons at 10~ 16 level will provide new tests in fundamental physics such as 
an improved measurement of Einstein's gravitational red-shift, a search for a 
possible anisotropy of the speed of light and a search for possible space-time 
variations of fundamental physical constants, similar to that described above 
in section 0] The current most precise measurement of the red-shift was 
made by the space mission Gravitational Probe A (GPA) with an accuracy 
of 7 x 10" 5 [55 . PHARAO/ACES will improve this test by a factor 30. By 
allowing worldwide comparison between distant clocks, operating with dif- 
ferent atomic species, ACES will play a major role in establishing new limits 
for variations of fundamental constants. 

Finally, PHARAO/ACES will be a pioneering cold atom experiment in 
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space. The PHARAO technology can be extended for the development of a 
new generation of high performance inertial sensors and clocks using mat- 
ter wave interferometry. As for atomic clocks, such sensors may achieve 
extremely high sensitivity in micro-gravity environment, as pointed out in 
the ESA HYPER project [56J. These instruments could then be used for a 
large variety of scientific space missions such as VLBI, gravitational wave 
detection, and deep space navigation. 

6 Conclusions 

With methods described in this paper, we expect to bring the accuracy of 
133 Cs fountains at 1 or 2 parts in 10 16 . For 87 Rb, a frequency stability of 

1 x 10~ 14 t -1 / 2 i.e. 3 x 10~ 17 at one day seems accessible, together with 
an excellent accuracy. Routine operation of these devices over several years 
will have a profound impact on ultra-precise time keeping and fundamental 
physics tests. To take full benefit of this performance, long distance time 
transfer systems must be upgraded. In particular, the ACES time and fre- 
quency transfer system will enable comparisons at the level of 10 -16 per day 
in 2008-2009. Another route currently under study makes use of telecom 
optical fibers and over 100 km distance a stability of 1 x 10 -14 at 1 s and 

2 x 10~ 17 at one 1 day has already been demonstrated jHZ]. Extension to 
larger distances is under study. 

More generally, clocks operating in the optical domain rather than in 
the microwave domain are making rapid progress on the ground [SB]- The 
frequency of these clocks is four to five orders of magnitude higher than the 
frequency of microwave standards and with an equivalent linewidth, the qual- 
ity factor of the resonance exceeds that of cesium clocks by the same factor. 
Using laser cooled atoms or ions and ultra-stable laser sources jSH], these 
optical clocks will likely open the 10~ 17 — 10~ 18 stability range. Using the 
wide frequency comb generated by femtosecond lasers, it is now possible to 
connect virtually all frequency standards together throughout the microwave 
to ultra-violet frequency domain jlHl EH] • The attractive proposal of [SH H>2| 
to realize an optical lattice clock is currently receiving a great deal of inter- 
est. In this method, neutral atoms are confined in an optical lattice in the 
Lamb-Dicke regime. Light-shifts of the clock levels induced by the lattice 
beams are differentially compensated at an appropriate laser detuning. This 
proposal combines several interesting features such as long observation time, 
large number of atoms, and recoil- free resonance [S3]. Promising atoms to 
implement this method are alkaline-earth atoms because of their strongly 
forbidden inter-combination line. Ca [IBUM], Sr [OBI EE] and Yb jBHJEIllHB] 
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are actively studied. 

In the frequency stability range of 1CT 17 — 1CT 18 , it is clear that fluc- 
tuations of the Earth potential at the clock location induced, for instance, 
by sea tides will affect comparisons between distant clocks. This limitation 
could be turned into an advantage if one installs such ultra-stable clock in 
space where the gravitational potential can present far reduced fluctuations 
compared to ground. As in the past, clocks with very high stability will have 
an ever increasing impact on scientific and industrial applications. 
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